Vetoing atmospheric neutrinos in a high energy neutrino telescope 
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Abstract 

We discuss the possibility to suppress downward atmospheric neutrinos in a high energy neutrino 
telescope. This can be achieved by vetoing the muon which is produced by the same parent meson 
decaying in the atmosphere. In principle, atmospheric neutrinos with energies E u > 10 TeV and 
zenith angle up to 60° can be vetoed with an efficiency of > 99%. Practical realization will depend 
on the depth of the neutrino telescope, on the muon veto efficiency and on the ability to identify 
downward moving neutrinos with a good energy estimation. 



Neutrino telescopes such as IceCube at the South Pole 111, ANTARES in the Mediter- 
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ranean Sea [2J and the Lake Baikal detector [3| search for neutrinos of extraterrestrial origin 
by using the Earth as a filter to suppress the background of atmospheric muons. The charged 
current channel is favored at TeV to PeV energies because of the large muon range, up 
to several kilometers of matter. To prevent downward-going atmospheric muons (j/i) from 
being mis-identified as muons induced by neutrino interactions, neutrino telescopes typically 
search for neutrino-induced upward or horizontal muons (t^ M -t/i) (Markov and Zheleznykh 

fl)- 

The principle of detecting neutrinos by looking for implies that the field of view of 
a neutrino telescope is limited to half of the sky i.e. the opposite hemisphere with respect 
to the geographical position of the detector. Moreover, with this approach, atmospheric 
neutrinos that penetrate through the entire Earth become an irreducible background for 
the search of extraterrestrial neutrinos. There are two standard approaches to separating 
extraterrestrial from atmospheric neutrinos in the sample. One is based on seeing an 
excess of events in a particular direction/time interval (point source search 5|), while the 
other is based on the assumption that astrophysical neutrinos have a harder spectrum than 
the atmospheric background (diffuse search {(J). 

The new generation of neutrino telescopes like IceCube [1], to be completed in 2011, 
and the R&D project KM3NeT will operate instrumented volumes of about one km 3 . 
Among the new opportunities offered by such large detectors, one feature seems especially 
interesting to us: the use of part of the instrumented volume as an active veto for l/j,. This 
opens the field of view of neutrino telescopes to the hemisphere above the detector jq]. 

At the energies involved in a neutrino telescope (TeV- PeV), the opening angle between 
the lu^ and the j/i produced in an atmospheric meson decay is very small. This implies 
that an atmospheric has a certain probability to arrive in the detector accompanied by 
its partner 

In this letter, we discuss the conditions under which a j/z-veto will also veto atmospheric 
ji^'s. This opens the opportunity to suppress what it is considered so far the irreducible 
background in neutrino telescopes. 

At high energy, neutrinos come predominantly from decay of charged kaons and pions. 
Charged pions decay with a probability of 99.99% through the reaction 7r ± — > /i ± + v^ipp) 
and charged kaons with a probability of 63.4 % via K ± — > ^ + v^iv^). There are small con- 



tributions from other channels, but we start by analyzing the simplest and most important 
case of pion and kaon decay and discuss the other cases at the end. 




FIG. 1: Two body decay of the parent meson into muon and neutrino. The left figure displays 
the back-to-back kinematics in the meson center of mass (cm) frame. The right figure shows the 
momenta after Lorentz transformation into the laboratory frame. 

Momentum conservation requires that the neutrino and companion muon are emitted 
back-to-back in the center of mass (cm) frame. The energies of the muon and neutrino in 
this frame are 
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(c = 1). Here T{ = m^/mf are the quadratic mass ratios, where the parent meson masses 
are m« (i = tt^, K ), the muon mass is and we neglect the neutrino mass. Numerically, 
r n = 0.573 and Tk = 0.046. Given the higher rest mass of the kaon with respect to that 
of the muon, the energy is quasi-equally shared between neutrino and muon for the kaon 
two-body decay, while for pion decay the energy balance is shifted more towards the muon. 
It should be noted that the muon energy is always larger than the neutrino energy in the 
cm-frame. 

Figure Q] displays the kinematical relations in the cm-frame and after Lorentz- 
transformation into the laboratory (lab) frame. After Lorentz-transformation along the 
positive x-axis from the cm- to the lab-frame, the muon and neutrino energies are given by: 

E v = 1 El m + (3 1P c ™ 

E» = 7E™ + f3 1P ™, (2) 
where 7 and (3 are the Lorentz factor and speed of the parent pion or kaon. 



Taking into account the back-to-back emission of neutrino and muon in the cm-frame, 



= \p cm \ cosO,, and (3) 
VZ = \p cm \ cos(0„ - tt) = -\p cm \ co&B v , 

where 6 U is the angle of the neutrino in the cm-frame of the parent meson relative to its 
direction in the lab-frame. In the approximation /3 — > 1, which is valid for meson energies 
above several GeV, we can then rewrite Eq. [2] as 

E v = 7|p cm |(l + cos^) and 

E, = 7b cm |(^-cos^). (4) 

The extreme condition for a given parent energy E{ = E v + E^ is obtained for cos 6 V = 1 
for which E^ min = riE U)max / (1 — r$). Thus, for any [E^E V ) the minimum energy of the 
companion muon is 

E, > (-^-) E v and 

Ei = E„ + E u > E v (—^—) (5) 

with % = tt , -ft' 1 ' 1 . The corresponding numerical values for pion and kaon decays are 

tt* : K > 2.342 ■ (6) 
K± : E K > 1.048 ■ ^ . (7) 

The corresponding numeric values for the muon-neutrino relations are E^ > 1.342 ■ E v 
for 7r ± decays, and E^ > 0.048 • E u for K . It should be noted that - unlike the situation 
in the cm-frame - the companion muon energy can be as low as 5% of the neutrino energy, 
if the parent meson is a kaon. The difference between pion and kaon decays which is related 
to the quadratic mass ratios r«, has an important impact on the energy dependence of the 
veto efficiency as discussed below. In particular, in addition to satisfying Eq. [51 the muon 
must have enough energy to penetrate to the detector. 

To veto [u^, the j/z track must also be relatively near to the neutrino trajectory. From the 
values of the transformed angles M , Q v we obtain typical distances between the neutrino 
and companion muon tracks after 10 km path length of less than 1 m (0.1 m) for neutrino 
energies above 1 TeV (10 TeV) if the parent meson was a pion, and less than 10 m (1 m) 



for kaons. High energy atmospheric neutrinos and their companion muons can therefore 
be treated as quasi-aligned given the typical granularity of optical modules in a neutrino 
telescope. 

The production spectrum of atmospheric neutrinos from tt^ (K) — > fi^ + is obtained 
from the convolution of the parent meson decay probability with the spectrum of mesons 
and the phase space distribution of the neutrinos: 



V V (E V ,X)= (8) 

U(E 1 X)dE 
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Here V u dE u is the number of neutrinos [y^ + v^) with energy between E v and E v + dE u 
produced per g/cm 2 along the direction defined by zenith angle 9. U(E,X) and K(E,X) 
are the differential energy spectrum of charged pions and kaons at slant depth X (in g/cm 2 ), 
and E is the energy of the parent pion or kaon. The factors in curly brackets are the decay 
probabilities of pions (kaons) at vertical depth Xcos(6). The factors in square brackets 
give the branching ratios times the normalized distributions of neutrino energies. The decay 
distributions are isotropic in the parent rest frame (uniform in cos^), so from Eq. HI the 
distributions of neutrino energy are flat over the kinematically allowed regions of phase 
space: < E v < E w (1 — tv) and < E v < E K (1 — r K ). The total differential intensity 
of neutrinos is obtained by integrating Eq. [8] over the whole atmosphere with E n ^ min and 
E K ,min given by Eqs. M and 

The contributions from pions and kaons in Eq. [H]are identical in form, but the kinematics 
are significantly different because of the difference in mass ratios. The critical energies below 
which decay is favored over hadronic re-interactions are also different, e n ~ 115 GeV and 
€k ~ 850 GeV. As a consequence of these differences the contribution of pions decreases 
and kaon decay becomes the main source of atmospheric neutrinos at high energy. (See n 
fraction in Fig. [2j ) 

If we require a minimum muon energy at production so that the muon can penetrate to 
the detector at slant-depth X, then, in addition to Eq. [5] 



E > E u + -E^minPO 



(9) 
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FIG. 2: Unaccompanied vertical atmospheric [v^ flux (upper); probability of accompaniment 
(lower) for seven depths. The grey curve is the fraction of [u^ from ir decay. 

must also be satisfied. The latter condition governs until E„ > E^ min (X) x (1 — r)/r, which 
is E v > 0.75 x E^miniX) f° r pi° n decay but E u > 20 x E^ min (X) for kaon decay. For 
sufficiently high neutrino energy, the accompanying muon is guaranteed at depth, but the 
asymptotic energy occurs later for the dominant kaon component. 

Fig. [2] illustrates the behavior of the proposed atmospheric veto as a function of depth 
at vertical incidence. The lower plot shows the probability that the partner muon reaches 
various depths as a function of E v . In each case the early onset of veto for the it — > \w 
component is indicated by the shoulder of the curve. The upper panel shows the remaining 
flux of atmospheric after the j/i veto is applied. Fig. [3] shows how the veto probability 
depends on zenith angle at two different depths, 1.8 and 3.5 km.w.e. The first approximates 
the center of IceCube at 2 km in ice, while the second represents the center of a detector at 
the NEMO site, which is the deepest candidate location for Km3NeT Ml. 

n 

For this illustration we use a simple energy-independent, average relation from MMC [9!] 
for the muon energy at production needed to reach depth X: 

E^winiX) = 0.73 TeV x {exp [X/2.8 km.w.e.] - 1 } . 

The veto probability decreases with depth for a given neutrino direction and with in- 
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FIG. 3: Veto probability as a function of zenith angle for 1800 and 3500 m.w.e. 



creasing zenith angle at a fixed vertical detector depth. These effects are both the result of 
the increased muon energy needed to penetrate to the deep detector. The dominant kaon 
contribution reaches its asymptotic value about a factor 30 higher in energy than the pion 
component. 

The values in Figs. [2] and [3] are obtained from an analytical approximation applicable for 
a power law primary cosmic-ray spectrum [lOj . For pions 
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with a corresponding expression for kaons. N (E) is the power-law primary spectrum of 
nucleons evaluated at the pion energy. Cross coupling between kaon and pion channels has 
been neglected as well as production of anti-nucleons. Z^n is the spectrum- weighted moment 
for pion production, A at is the nucleon interaction length and Ajv, A„-, are attenuation 
lengths for nucleons, charged pions and charged kaons, respectively. 

Eq. M gives the production spectrum of neutrinos as a function of slant depth X in the at- 
mosphere. The total differential neutrino spectrum requires evaluating j§ round V U (E, X)dX. 
High-energy muons originate high in the atmosphere, so it is a good approximation to ex- 



tend the upper limit of this integral to infinity. The result for the contribution from charged 
pions is 

MEu) = N {E v )Z N ^_m x joo. n _^ (n) 

1 Att/Ajv-I i (A n /A N -1) 2 , 

_z+£{E v ) 2z+(,(E v ) 3z+£(E„) I" • • -J > 

with a similar expression for kaons. Here z = E / E v and z m i n is the greater of 1 + -E M , m in / E v 
from Eq. [9]or 1/(1 - r<) (Eq. [5]). Also, = e t /Ecos(9). 

The integral in Eq. [TT]can be evaluated analytically in the limits of low {E v cos(#) <C e n ) 
and high (i±/V cos(#) ^> e n ) energy. One can then combine the low and high energy limits 
into a single approximation, 
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The expression [12] with numerical values of the parameters from Ref. [10|] is used to make 
Figs. Hand El 

The new opportunity to veto high energy atmospheric muons using part of a neutrino tele- 
scope js], opens the possibility to suppress downward going atmospheric neutrinos. Extra- 
terrestrial neutrinos, which are in general not accompanied by a muon, will not be discarded 
by the veto system. If a downward event starts inside a fiducial volume of the detector and 
has no activity in the outer veto region, and has sufficiently high energy that an atmospheric 
[v^ would be accompanied by its partner muon, then the neutrino can be classed (with a 
certain probability) as being of extra-terrestrial origin. 

So far we have assumed that all neutrinos come from decay of charged kaons and pi- 
ons. We now estimate the contribution of atmospheric ji^ from minor channels and assess 
their effects on the veto probability. The semi-leptonic decay has a branching ratio of 
0.27 (K®) and 0.033 {K^). The final state is n, //, so the minimal kinematic configuration 
occurs as before when the muon is backward in the cm system of the decaying kaon, balanced 
by the forward moving neutrino and pion. In this case, however, the forward momentum 
is shared by two particles so a lower E v guarantees the presence of the partner muon for a 
[Py,. This channel therefore does not weaken the veto probability. At the other extreme are 



the from muon decay, in which case there is no muon. The contribution is small at high 



For E v ~ 10 TeV and zenith angles < 60°, this 



energy and can be estimated as in Ref. 
contribution is less than one per mil. 

Contributions of prompt neutrinos from charm are more complex and more difficult to 
assess. A model for charm production that nearly saturates existing limits 12| is the Re- 
combination Quark Parton Model (RQPM) [13]. The charm contribution to atmospheric 
\.Vfj, is less than 10% for E v = 10 TeV in this model and crosses over the conventional 
contribution just above 100 TeV. Relevant decay channels include D + — > K yU + (7%), 
D + — > K~ 7r + fi + (4%), and A+ — > Afi + (2%). Because the charmed mass is higher, the 
energy fraction carried by the muon can be even lower than in the case of the kaon. On the 
other hand, the final states involve several particles so the neutrino energy required to see 
the partner muon is relatively low as well. Given the branching ratios and the upper limits 
on charm contribution, unaccompanied prompt should weaken the veto probability by 
less than one percent. 

We conclude therefore that it should be possible in principle to veto atmospheric neutrinos 
with energies in the multi-TeV range and zenith angles less than 60° with an efficiency of 99% 
or somewhat better. This paves the way for a sensitive and new type of search of neutrinos 
from astronomical sources. The main limitation will be the extent to which real detectors can 
put a lower limit on the energy of a neutrino interaction that starts in the detector. Because 
of the steep spectrum of atmospheric neutrinos, most of the interactions will be near the 
nominal threshold, thus enhancing the smearing effect of fluctuations. Detailed simulations 
of individual detectors, including surface air shower arrays, or alternative veto systems at 
shallow depths, are needed to assess the veto capability in practice. Such simulations will 
also be able to include the case where the [v^ is vetoed by a muon on a different branch in 
the accompanying air shower, which enhances the veto probability to some extent at high 
energy. 
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